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We analyse recently proposed physical implementations of a quantum computer based on polar 
molecules. A set of general requirements for a molecular system is presented, which would provide 
an optimal combination of quantum gate times, coherence times, number of operations, high gate 
accuracy and experimental feasibility. We proceed with a detailed analysis of a scheme utilizing 
switchable dipole-dipole interactions between polar molecules. Switchable dipole-dipole interaction 
is an efficient tool for realization of two-qubit quantum gates, necessary to construct a universal set 
of gates. We consider three possible realizations of a phase gate using specific molecules, such as 
CO, NF, alkali dimers and alkaline-earth monohalides. We suggest suitable electronic states and 
ransitions and investigate requirements for the pulses driving them. Finally, we analyse possible 
sources of decoherence. 

PACS numbers: 



I. INTRODUCTION 

Information processing relying on quantum mechani- 
cal rather than classical systems holds the promise for 
a dramatic speedup of such operations as factoring large 
numbers, searches of unstructured databases and simula- 
tion of the dynamics of quantum mechanical systems [l| . 
Communication using optical fields with single photons 
would provide a level of security impossible with clas- 
sical communication techniques (2j. Recent years have 
witnessed remarkable advances in both the theoretical 
and experimental development of quantum computing 
technologies, including demonstrations of basic building 
blocks necessary for quantum computing and quantum 
networking. Various approaches have been explored, but 
the most advanced are based on trapped ions and neu- 
tral atoms d-d, cavity QED [|, liquid NMR @|, and 
solid-state systems Q- Recently dipolar molecules were 
proposed as a system with a set of parameters optimal for 
physical implementation of quantum computing schemes 
Q. Polar molecules combine the advantages of neutral 
atoms and ions (such as long coherence times, rich level 
structure, strong optical and microwave transitions, well- 
developed techniques of coherent manipulation with op- 
tical and microwave pulses) and of quantum dots and 
superconducting circuits (easy control with electrostatic 
fields). Polar molecules are thus compatible with var- 
ious architectures, including optical lattices, microwave 
and electrostatic traps and solid-state systems. An im- 
portant aspect of these systems is the ability to access 
electronic states which exibit a large permanent dipole 
moment. The molecules can then be used for fast con- 
ditional dipole-dipole interactions resulting in two-qubit 
operations, necessary for construction of a universal set 
of qubit gates. 

A number of specific implementations of quantum com- 
puting schemes with polar molecules have been sug- 
gested. In the original proposal [8| the two projections of 



a permanent dipole moment of a polar molecule on the 
direction of an external electric field are used to store a 
qubit. Molecules are held in a ID trap array with a gradi- 
ent of the electric field, producing a Stark shift specific to 
each molecule, thus allowing molecules to be addressed 
individually by spectroscopic means. A two-qubit gate 
is implemented via electric dipole-dipole interaction of 
molecules. In Ref. [9| a qubit is encoded in a bound molec- 
ular and a free state of two atoms in an optical lattice site. 
Free atoms can be transferred into the bound molecular 
state and back with a Raman pulse. If the bound molecu- 
lar state has a large permanent dipole moment, then two 
molecules in neighboring sites will dipole-dipole interact 
producing a phase shift. The phase shift accumulated by 
a molecule is conditional on the state of another molecule, 
thus resulting in the phase gate. Ref. [HI further develops 
the idea of Ref.[8[ for polar molecules in optical lattices. 
Tunable dipole-dipole interaction between two molecules 
can be realized when a microwave field resonant with the 
J = — > J = 1 rotational transition in the ground elec- 
tronic and vibrational state induces a dipole moment in 
each molecule. 

Based on this, a set of general requirements for a dipo- 
lar molecular system is proposed, which would provide 
an optimal combination of long coherence times, short 
gate times (resulting in maximal number of operations), 
small gate errors, and experimental feasibility. 

1) Choice of qubit states 

To store a qubit long-lived states are required, well 
isolated from the environment, i.e. weakly perturbed by 
electric and magnetic fields and various interactions (for 
example, dipole-dipole or spin-spin). Good candidates 
are hyperfine and rotational states of a ground electronic 
molecular state having a negligible permanent dipole mo- 
ment. 

2) Coupling strengths 

Fast one and two-qubit gates require the correspond- 
ing interaction strengths to be large. The storage states 
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need long lifetimes. This means that the transition be- 
tween the qubit states |0) and |1) is forbidden and thus 
Raman transitions could be used to perform one-qubit 
gates. An alternative approach is to map the qubit to 
some states coupled by a one-photon allowed transition. 
In both cases, the resulting transition coupling has to be 
strong. Choosing electric dipole-allowed transitions with 
dipole moments of a fraction of a Debye for one-qubit 
operations, and sufficiently large Rabi frequencies of the 
laser fields, one qubit gate times can be as small as fem- 
tosecond laser pulse durations. 

3) Robust dipole-dipolc interactions 

Possible errors resulting from dipole-diplc interactions 
include (i) instabilities in inter-molecular distance and 
alignment, and (ii) instabilities of the dipole-moment 
states themselves. To avoid (i), we need either "frozen" 
molecules or a more robust dipole-dipole coupling mech- 
anism such as dipole blockade Permanent dipole 
states with life-times long compared to the gate times 
are necessary against (ii). 

4) Cooling and trapping 

In order to fulfill the requirements of points (2) and (3), 
the molecule has to be cooled down to sub-Kelvin tem- 
peratures (typical rotational transition frequencies ~ 10 
GHz, giving the condition T < 1 K to avoid population 
of higher-energy rotational states) . In optical lattices and 
electrostatic traps it is necessary to have the molecules 
occupy a ground motional state of the trap to minimize 
decoherence and gate error. 

5) Decoherence 

Provided that a qubit is encoded in hyperfine or ro- 
tational sublevels with typical lifetimes of the order of 
hours and typical optical lattice or other trap decoher- 
ence times on the order of seconds, the main decoherence 
mechanisms are (i) limited lifetimes of the (excited) high- 
dipole states, (ii) spatial dependence of the dipole-dipole 
interaction, and (iii) finite laser-linewidth and other light- 
induced decoherence. Points (i) and (ii) are discussed 
above, and we expect decoherence times stemming from 
finite laser linewidth of the order of ms using phase- 
stabilized lasers [121 ]. 

Recently we proposed a model of switchable dipole- 
dipole interaction between polar molecules [l3j], allow- 
ing one to implement a universal two-qubit gate, such 
as a phase gate. In the proposed scheme, the state used 
to store a qubit and an additional state to switch the 
interaction have significantly different permanent dipole 
moments. Ideally, the storage state has a zero moment, 
while the " switch" one has a large dipole moment of sev- 
eral Debye, so that two molecules excited to this state 
will exibit strong dipole-dipole interaction and acquire a 
7r phase shift. Provided that the molecules can be ex- 
cited only from one of the qubit states, for example |1), 
only the |11) state of the two-qubit system flips sign, 
resulting in a phase gate. In this approach only two se- 
lected molecules interact, which drastically simplifies the 
two-qubit gate compared to the case when all molecules 
interact at the same time. In this work we carry out a de- 



tailed study of the scheme considering specific molecules 
as examples, and check if it satisfies the general require- 
ments presented above. 

Three realizations of the phase gate with polar 
molecules were proposed in previous work. The first one, 
a direct phase gate, uses a molecular state with a zero 
permanent dipole moment to store a qubit, while to per- 
form the gate operation, molecules are excited to an ad- 
ditional state with a large dipole moment (see Fig. (QJ,)). 
We also consider the "inverted" case shown in Fig. Hb), 
in which the storage state has a large dipole moment, and 
the additional state e) has a small one. In this case the 
molecules interact and acquire a phase shift when both 
are in the storage state. The two systems can be im- 
plemented using different electronic states of a molecule, 
as will be shown below. In a "rotational" scheme, we 
could realize both the direct and inverted schemes using 
rotational levels of the ground state rather than excited 
electronic states. The direct and inverted schemes using 
electronic molecular transitions in the visible or UV range 
can be most naturally implemented with cold molecules 
in optical lattices [3[. The direct case can also be real- 
ized with molecules doped into solid-state matrices. The 
rotational scheme can be implemented using a recently 
proposed architecture combining electrostatic traps and 
coupled microwave resonators [14|. 

The organization of the remainder of the paper is as 
follows. In Section II we analyse the constraints imposed 
on the relative distance and angles of the two interacting 
molecular dipoles by the requirement to keep a phase er- 
ror below a threshold value (we set it at 1%). Section III 
considers the direct dipole switching scheme with a CO 
molecule. The inverted scheme is analysed in Section IV. 
The scheme based on rotational states is studied in Sec- 
tion V. Decoherence mechanisms are analysed in Section 
VI. Finally, we conclude in Section VII. 

II. PHASE GATE ERROR DUE TO 
MOLECULAR MOVEMENT AND DIPOLE 
MISALIGNMENT 

We analyse the phase errors assuming that two inter- 
acting molecules are transferred into states \e) having 
large permanent dipole moments, and the phase is ac- 
cumulated while they undergo dipole-dipole interaction. 
The phase dependence on the relative distance and ori- 
entation of the dipole moments is given by the expression 

, _ T ( 3Qhf)(jhf) _ fh_jh\ = 
Ti \ r 5 r 3 J 

T u 2 ( 

= =- 1 3 sm(0 + 9\ ) (cos #2 sin + sin O2 cos 9 sin (f>2 ) 

h r 6 \ 

— cos 6i cos 02 — sin Q\ sin 62 sin , (1) 

where the angles 9, 9\, #2 and (f>2 determining the orien- 
tations of the dipole moments fl± and /I2 of two interact- 
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FIG. 1: Schematic showing direct (a) and inverted (b) vari- 
ants of the phase gate based on controllable dipole-dipole in- 
teraction of polar molecules. 



ing molecules are shown in Fig[2J T is the duration of the 
gate operation, r is the vector connecting the first dipole 
with the second one. The vector of the first dipole jli is 
assumed to be in the (y,z) plane as well as the r vector to 
simplify the analysis (we can always choose the coordi- 
nate system this way). The angle 9 describes the vertical 
shift of the second dipole with respect to the first one. 
The angles &i and 02 are the polar angles of the dipoles 
in their coordinate systems, 4>2 is the azimuthal angle of 
the second dipole. We will analyse the sensitivity of the 
accumulated phase to the relative distance between the 
molecules and the angles, assuming small deviations from 
their equilibrium values r = (r), = 0\ = 02 = 0, and 

4>2 = 7T/2. 

Fidelity of the phase gate, which is directly related 
to a phase error during the gate operation, will then be 
affected by the relative motion of the molecules, mis- 
alignment of their dipole moments, and stability of the 
duration of optical excitation pulses. From Eq.([T|) the 
phase error can be expressed via the uncertainty in the 
distance r and angles 0, 0\, 02 and 02 ■ 
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Requiring the phase error be < 1% we see from Eq.Q 
th at the u ncertainty of the relative distance has to be 
^/((Ar) 2 ) < 0.003 (r). That is for (r) - 500 nm the 
molecules have to stay within 1.5 nm of their average 
distance. As follows from Eqs.([3])-((5]) the same phase 
error of 1% results in V((A6>) 4 ) ~ 3 • 10~ 3 (or A0 - 3°), 
and ^((Atf^) 4 } ~ 2 ■ 10~ 2 (or A0 1)2 ~ 8°). 

For a simple estimate of the mean distance a molecule 
travels in an optical lattice potential minimum we make 
a harmonic approximation of the potential Vq sin 2 (fca;) w 
Vok 2 x 2 , giving the frequency uj = k^2Va/m with the cor- 
responding motional ground state wave function width 
a = yjh/mu). Taking V = (10 - 40)E R = (10 - 
A0)h 2 k 2 /2m, where Er = h 2 k 2 /2m is the recoil energy 
and (10 — 40)E R is the typical depth of a lattice potential, 
the width a can be related to the lattice field wavelength 
A and the lattice depth as a = (10 - 40)" 1/4 A/27r. For 
molecules in neighboring lattice sites (r) = A/2, and the 
phase error according to Eq.([2]) is 6(10 — 40) _1 / 4 /27r rs 
0.4 — 0.55, for five lattice periods distance the phase error 
is reduced to ~ 10%. To reduce the error even further the 
mean separation between the molecules in the lattice can 
be increased along with the depth of the potential. At 
the same time this will lead to a decrease in the strength 
of the dipole-dipole interaction and, therefore, to a larger 
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T of the phase gate. 

A better solution is to make use of the dipole blockade 
phenomenon [llj]. In this case the excitation scheme B 
of Ref.[ll| can be applied, in which the doubly excited 
state |ee) of the two molecules is shifted due to the dipole- 
dipole interaction from its unperturbed position, the shift 
being u ~ fj? /r 3 h, where fi is the permanent dipole mo- 
ment of the \e) state. To implement the diplo-blockade 
mechanism the two molecules have to be addressed in- 
dividually. Then, the first part of the gate operation is 
the transfer of a control molecule to the excited state |e) 
with a resonant tt pulse. The second step is to apply to 
a target molecule a 2tt pulse, resonant with the unper- 
turbed transition |1) — -» |e). Since the |ee) state is de- 
tuned, the target molecule returns back to |1) acquiring 
a phase 4> ~ t^^-k/^u <C 7r if u f^- In the last step 
the control molecule is brought back to 1 1) with a second 
7r pulse. The major requirement of the dipole-blockade 
mechanism is to keep the Rabi frequency of the 2tt pulse 
much smaller than the energy shift of the doubly excited 
state. As a result all states except the |00) one acquire a 
7r phase shift, resulting in the phase gate. 

Since the exact position of the two molecules is not 
important if the dipole blockade mechanism is used, it 
allows one to circumvent the strict requirement on the 
relative distance set by Eq.([2]) and on angles by Eqs.Q- 
^ . The smaller the uncertainties in the relative distance 
and angles, the better the dipole blockade mechanism 
works, resulting in smaller fluctuations of u. The an- 
gle 9 can be made small by loading molecules in a 3D 
lattice, which would provide confinement in the vertical 
direction. Combination of the vertical confinement with 
the dipole-blockade mechanism allows one to satisfy the 
9 angle requirement. The 9\^ angles can be controlled 
with either static or microwave electric fields. Finally, 
the (f>2 dependence as is seen from Eq. (O of the phase is 
weak, provided that the other two angles are close to the 
optimal values. 

The same analysis applies to the architecture utilizing 
electrostatic traps combined with coupled microwave res- 
onators [l4[ , where the dipole-dipole interaction strength 
scales as /i 2 /h 2 r with h and r being the distance between 
the molecule and the trap surface and the distance be- 
tween the traps, respectively. In this case the phase error 
dependence on the geometry of the trap is given by 



W (h) ' 

(<f>) (r) ' 

For a typical value of h ~ 1 /im the uncertainty of h 
has to be \J ((Ah) 2 ) < 5 nm to have the phase error be- 
low 1%, and with typical (r) ~ 10 //m the corresponding 
uncertainty */ ((Ar) 2 ) < 100 nm. Again, the require- 
ment for h is hard to fulfill unless the dipole blockade 
mechanism is used. 



In the next three sections we will describe three imple- 
mentations of the phase gate using specific molecules as 
examples. 



III. DIRECT PHASE GATE ON THE BASIS OF 
CO MOLECULE 

CO molecule has a small permanent dipole moment 
(0.1 D) in the ground electronic state X which 
makes it suitable for the direct phase gate implemen- 
tation. We choose as the \e) state the metastable a 3 ±Io 
electronic state with a permanent dipole moment of 1.37 
D 15]. If we choose an isotopic variant of CO, the qubit 
can be encoded in hyperfine sublevels of the ground rovi- 
brational state, weakly coupled with other states in the 
presence of static and non-resonant dynamic electromag- 
netic fields. We start with an analysis of the hyperfine 
structure of isotopic CO molecules in the X 1 S + and 
a 3 n states. 



A. Hyperfine structure of CO 

There are three CO isotopomers 13 CO (1% abun- 
dance), C 17 (0.038% abundance), 13 C 17 (3.8 • 10~ 4 % 
abundance) with at least one nucleus having a non- 
zero spin. The 13 C and 17 nuclear spins are 7c — 
1/2 and Io = 5/2, respectively. For the ground 
rotational state J = of X 1 E + electronic state 
the coupling with the nuclear spin is described by 
the Hamiltonian Hhf s = bJI — eQqF(I, J), where 
F(I,J) = (3C(C + l)/4- 7(7 + 1) J (J + 1)) /27(27 - 
1)(2J - 1)(2J + 3) is the Casimir function (with C = 
F(F+l)-I(I+l)-J(J+l)). For J = (F = I) the cou- 
pling term vanishes, resulting in zero hyperfine splitting 
in the ground state for both 13 CO and C 17 0. Small hy- 
perfine splitting is present in 13 C 17 due to interaction of 
the carbon and oxygen nuclear spins. The level structure 
of the lowest rotational states of the X 1 S + and a 3 Hq 
states is shown in FigJ3] [lfj. In the case of 13 CO the 
±1/2 nuclear spin sublevels of the ground state, Zeeman 
split by an external magnetic field, can be used as |0), |1) 
qubit states, respectively. As the e) state we can choose 
a Zeeman sublevel of J = 1, F = 1/2 hyperfine level of 
the long-lived a 3 IIo (v = 0, J — 1) state with a lifetime 
~ 0.5 s. Zeeman splitting in the ground electronic state 
is small, scaling as ~ 1 kHz/G, while in the excited state 
it scales as ~ 1 MHz/G due to the non-zero electronic 
angular momentum of the a 3 II state [lTJ • It means that 
magnetic fields ~ 10 G will suffice to make the Zeeman 
levels of the |e) state resolvable. Selective excitation of 
|1) can then be realized with a a + laser pulse resonant 
with the —1/2 — > +1/2 transition between the ground 
and excited electronic states. Single qubit rotations can 
be performed with Raman pulses using the same +1/2 
sublevel of the \e) state. 
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In the case of 13 C O the hyperfine sublevels F\ = 
2, 3 of the ground rovibrational state can be utilized as 
the |0), |1) states, respectively. As the |e) state we can 
choose the F\ = 2 component of the J — 1,F = 5/2 
hyperfine level of the excited state. Hyperfine splittings 
of ~ 10 — 100 MHz are expected in the a 3 U state due 
to strong electron spin - nuclear spin interaction, so that 
the hyperfine structure in this state is expected to be 
resolved with a narrow-band laser. Selective excitation 
of 1 1) can be realized then using a cr~ polarized laser pulse 
resonant with the J = 0, F = 5/2, Fi=3->J=1,F = 
5/2, F± = 2 transition, as shown in FigJ31 Single qubit 
manipulation can also be performed via this excited state 
hyperfine sublevel using a a~ and linearly polarized laser 
pulses. 

Readout of the qubit states can be done in the same 
manner using the short-lived (9 ns lifetime) excited elec- 
tronic A 1 ni(J = 1,F = 1/2) state, with the corre- 
sponding transition wavelength 147 nm. It is the first 
excited singlet state and therefore it decays directly to 
the X 1 S + state. Initialization of the qubit into the |0) 
state can be done by optical pumping using a CW a + (in 
the case of 13 CO) or a~ (in the case of 13 C 17 0) polarized 
laser beam applied in the same way as in Fig|3] 
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B. Phase gate operation time 



The spin-forbidden X 1 E+ 



3 n transition is weakly 



allowed due to the mixing of the A and a 3 Hi states 
by spin-orbit interaction and rotational mixing between 
a 3 Ui and a 3 n states [Uj]. The lifetime of the J = 1 
rotational level of the a 3 n state is 0.5 s [lj|, giving the 
effective transition dipole moment fii n d ~ 2 • 10~ 4 D. If 
we require that the n pulse, transferring population to 
the a 3 n state has duration T n — %/Q, ~ 50 /j,s, the cor- 
responding Rabi frequency and electric field amplitude of 
the pulse are Q.^ = 6 • 10 4 s -1 and E ~ 10 2 V/cm, with 
the intensity of the laser pulse I = cE 2 /4ir ~ 25 W/cm 2 . 

As was shown in the previous section the dipole block- 
ade mechanism has to be used to keep the phase error 
small. For two molecules in neighboring optical lattice 
sites with r = A/2 = 100 nm (assuming that molecules 
are trapped by a blue-detuned lattice field near resonant 
with the X 1 E + — a 3 n transition) the shift of the |ee) 
state due to dipole-dipole interaction is u = 1.87 • 10 6 
s _1 (« 300 kHz), so the gate can be performed with two 
7r pulses 7r c exciting and de-exciting a control molecule, 
and a 27r pulse 27r 4 applied to a target molecule. Choos- 
ing a Rabi frequency ~ 10 5 s _1 for both the it and 
2ir pulses (satisfying <c u) results in the gate time 
Tgate — In/Sin + 2h/Q,2t; ~ 126 /is. If the molecules 
are separated by five lattice periods, for example, the en- 
ergy shift of the |ee) state is u rs 1.5 • 10 4 s _1 (« 2.4 
kHz), the Rabi frequency has to be reduced to fulfill 

2 • 10 3 s _1 , the gate time 
« 6.3 ms. Since the dipole- 
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FIG. 3: Hyperfine structure and Zeeman splittings of the low- 
est rotational states of the ground X J E + and excited a 3 Lto 
electronic states of two isotopomers of CO along with schemes 
of selective excitation of |1) — > |e) transition. 
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sitions of the molecules is not affecting the gate. The 
resulting phase gate has the following form: 
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Another candidate for the direct phase gate is the NF 
molecule, which has a small dipole moment fi w 0.075 D 
in the ground electronic X 3 £~ state and \i « 0.75 D in 
the excited metastable b 1 E + state (2lj . 



IV. INVERTED PHASE GATE WITH LiCs 

The inverted phase gate can be implemented in mixed 
alkali dimers, many of which have a large dipole moment 
in their ground electronic state X 1 E + and a small dipole 
moment in their metastable a 3 E + state. Heteronuclear 
alkali dimers are of growing interest due to studies of 
cold collision dynamics and photoassociation processes 
of laser cooled alkali atoms. We analyze the inverted 
gate usin g Li Cs, which has recently been created experi- 
mentally [25(| from the ultracold atomic gases Li and Cs. 



A. Hyperflne structure of LiCs 

The electronic structure of LiCs was calculated in [2(| , 
and a schematic figure of the six lowest electronic states 
is approximately drawn in Fig0J The metastable a 3 £ + 
state supports a number of bound states (up to 20 vibra- 
tional states were observed in NaCs 27]). The ground 
X 1 S + state of LiCs has a dipole moment of fss 5.5 D and 
the metastable state has a permanent dipole moment of 
w -0.45 D H|. 

Taking into account that the 7 Li has a nuclear spin 
In = 3/2, and 133 Cs has a nuclear spin Ic s = 7/2 
the projected ground state X 1 Y>{v = 0,J = 0) hyper- 
fine structure along with the experimentally observed 
metastable state a 3 E + hyperfine structure, are shown 
schematically in Figj5] The ground state hyperfine struc- 
ture has not yet been observed in mixed alkali molecules, 
and the splittings are most probably in the kHz - hun- 
dreds kHz range. The hyperfine structure of the c 3 S + 
and a 3 S + states was observed in 23 Na 85 Rb and in 
23 Na 39 K with typical splittings of F2 and G2 levels 
- 30 - 40 MHz and - 300 MHz [M H H, respec- 
tively. For LiCs we assume splittings of the same order. 
Hyperfine sublevels F — 4 and F = 5 of the ground 
rovibrational state can be chosen to encode the qubit |0) 
and |1) states, respectively. The excited |e) state can be 
encoded in the G± = 9/2, G2 = 6 hyperfine sublevel of 
the metastable a 3 S + (u = 0, N = 0, J = 1) state. Using 
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FIG. 4: Schematic of the low-lying electronic states of LiCs. 



circularly polarized laser pulses the |1) state can be se- 
lectively excited into |e) as shown in FigfSb, since the G2 
sublevels are well resolved. 

The spin- forbidden X 1 S + — > a 3 E + transition is 
weakly electric-dipole allowed due to spin-orbit mixing 
of the a 3 S + state with high energy singlet electronic 
states. In order to obtain an estimate of an effective 
dipole moment of the transition we can use as exam- 
ples the known parameters of spin-forbidden transitions 
in other molecules. Weak fluorescence of the a 3 S - 
X 1 E + was observed in BiN, BiP and BiAs [36| with ra- 
diative lifetimes — 0.5 ms. The CO molecule, considered 
in the previous section, has the spin- forbidden X 1 E + - 
a 3 n transition with lifetimes ms - hundreds ms depend- 
ing on the fine and rotational level. If we assume the 
same lifetime of a 3 E + for LiCs as for BiN, and take into 
account the wavelength of the transition ~ 2 /xm, the 
corresponding effective transition electric dipole moment 
is Hind ~ 0.2 D. 



B. Phase gate operation time 

Cold alkali dimers are currently produced from cold 
alkali atoms usin g ph otoassociation and Feshbach reso- 
nance techniques [23, H3, HH]. Optical lattices with one 
molecule per site have been demonstrated [30l |. Thus 
alkali dimers are natural candidates for quantum com- 
puting in optical lattices. As was shown above the phase 
accumulated during the phase gate strongly depends on 
the positions of the interacting molecules, and it is nec- 
essary to utilize the dipole blockade mechanism. The 
inverted gate can be slightly modified to use it. Let us 
assume that the molecules can be addressed individually. 
Then in the first step we excite the control molecule from 
X = 0, J = 0,F = 5) (|1) state) to some higher- 

energy vibrational state X 1 E + (w = vq, J — 0, F = 5) 
(|1') state), having significant overlap with the ground 
vibrational state of the a 3 E + (w = 0, N = Q, J = 1) 
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FIG. 5: Hyperfine structure in the ground X 1 E + (w = 0, J = 
0), X 1 E+(« = Wo, J = 0) and excited a 3 E+(t> = 0,iV = 
0, J = 1), c 3 E+(« = vi,JV = J = 1) states. Hyperfine 
coupling for c 3 E + and a 3 E + was found to be Hund's bpj 
and bps, respectively [32l ]. a) Population transfer from |1) 
state (ti = 0, J = 0, F = 5) to a higher vibrational level |1)' 
(v — vo,J = 0, F = 5) of the ground state, having optimal 
Franck-Condon factor with a ground rovibrational state of 
a 3 E + ; b) Selective excitation of |1)' into |e), encoded into 
a 3 Z+(v = 0, JV = 0, J = 1, Gi = 9/2, G 2 = 6). 



(|e) state). It can be realized using a STIRAP technique 
in the way shown in Fig(S^i. A it pulse next transfers 
the control molecule from to the lowest rovibrational 
state of the a 3 S + state (|e) state). 

Next the target molecule is Raman-excited into |1'), 
and a 2ir pulse resonant to the transition |1') — > |e) is 
applied. The control molecule is then deexcited with a 
second ir pulse into |1'), and finally both molecules are 
Raman transferred back into |1). If the control molecule 



is in the |0) state, it is not excited to |e), and interacts 
strongly with the target molecule, shifting the |1') — » 
\e) transition frequency. The probability of excitation 
of the target molecule to the |e) state during the 2tt is 
therefore very small. It will result in a negligible phase 
shift, significantly less than n, accumulated during the 
2n pulse, thus implementing the phase gate. 

The total gate operation time comes from the STI- 
RAP sequence and the w pulse for the control molecule, 
followed by the STIRAP sequence and the 2tt pulse for 
the target molecule and the STIRAP deexcitation. Fi- 
nal contribution is from the control molecule tt pulse 
and STIRAP deexcitation sequence. The intermediate 
c 3 S+ state lifetime is - 20 - 27 ns, so the STIRAP 
pulses have to be of 100 ps-1 ns duration. The dura- 
tion of the 7T pulse between the X 1 £ + (t>o, J = 0) and 
a 3 £+0 = 0, N = 0, J = 1) states can be set at T x = 50 
fj,s as for the CO case, with the 27r pulse duration of 100 
/is, respectively. For the 50 fis tt pulse the required Rabi 
frequency is f2 ff = 6 • 10 4 s _1 , resulting in the amplitude 
of the electric component of the laser pulse E ~ 0.15 
V/cm, with a corresponding laser pulse intensity / ~ 70 
/iW/cm 2 . The 27r pulse Rabi frequency has to satisfy 
two conditions: first, C u to be able to use the 
dipole blockade mechanism due to the interaction of the 
molecules in the ground electronic state; second, it should 
take into account that the dipole moment of the |e) state 
in the case of LiCs is not really zero. The Rabi fre- 
quency has to exceed the shift of the | el) state due to 
dipole-dipole interaction between the control molecule in 
\e) and the target molecule in the |1) state to have the 
27T pulse resonant. Taking the \e) state dipole moment 
w 0.45 D, and assuming that a blue-detuned lattice laser 
field for LiCs will have a wavelength A ~ 800 nm, for 
two molecules in neighboring sites the energy shift of the 
|1) — * |e) transition due to dipole-dipole interaction is 
M e i ~ 4 • 10 4 s" 1 and the 2-7T Rabi frequency of I^tt = 10 5 
s _1 will suffice. On the other hand the energy shift of the 
|1) — > |e) transition due to dipole-dipole interaction of the 
molecules in the ground electronic state is u ~ 5 • 10 5 s _1 
and u 3> r^w, allowing one to use the dipole-blockade 
mechanism. The STIRAP pulses are much shorter than 
the 7T and 27r pulses, so their contribution to the phase 
gate is small, and the total phase gate operation time is 

Tgate = 271-/^ + 2k/Q, 2 it ~ 160 /iS. 

From the rotational constant of the ground electronic 
state of LiCs B = 0.1935 cm -1 and its dipole moment 
/j, [26| the magnitude of a static electric field necessary 
to efficiently mix the J = and J = 1 rotational states 
of the ground state can be calculated £ ~ 2Bh//i w 4 
kV/cm. The electric field is switched on only during the 
27T pulse. 

The gate we obtain after the sequence of pulses de- 
scribed above is: 
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110} 


2ir t 


e i(*e + * t ) | 10 ) ffj 


e »(*c+* t ) | 10 ) 5 


111) ^ 


i|le> 


2tt ( 


-ie***/ 2 |le) ^ 


e l *'/ 2 |ll), 



(8) 

where in the third line of Eq.® we made use of the 
detuning of the state 1 10) from its unperturbed position 
due to the dipole-dipole interaction, so that the 2n pulse 
applied to the target molecule is not resonant and does 
not change its phase. In the last line of Eq.© line we 
took into account that during a resonant 2tt pulse the 
target molecule spends half of the pulse duration in |1), 
and therefore, accumulates a phase shift $t/2. The phase 
shifts $t,c and <E>t are accumulated during the 2ir pulse 
due to dipole-dipole interaction with all molecules in the 
lattice and due to evolution in the DC electric field: 



E 



fxST 2w /h, 



(9) 



jjtc,t C >*J 

fi 2 T y?T 



E 



[lETjh + Me 5T 27r /fi(10) 



where fx and fi e are the X 1 £ + and a 3 E + states per- 
manent electric dipole moments. 

Choosing £ sufficiently large so that the electric field 
induced phase is much greater than the interaction phase, 
we can set $ c ,t s» $t = ir(2n + 1), resulting in the gate: 



5/2 



-5 201.6 MHz 
-4 101 MHz 



3.2 GHz 



3/2 



xV 



3/2 



1.6 GHz 



1/2 




N 



rr 

I, >_ 
-V — 



-20.5 MHz 



-1 -34.3 MHz 

"2 -40 MHz 

-2 -64.6 MHz 

-1 -80 MHz 

-0 -100 MHz 

" 3 -110 MHz 



i — 4 -170.6 MHz 
i—i 4 167.4 MHz 

n 



-3 29.3 MHz 



-2 -52 MHz 
-2 -74.3 MHz 



-3-93 MHz 
-1 -143.3 MHz 



1/2 
.1 



.3 138 MHz 
.2 -193 MHz 



FIG. 6: Proposed scheme of selective excitation of |1 > state 
(encoded in the F=3 hyperfine sublevel of the N=0 rotational 
level of the ground vibrational state of X 2 S + electronic level 
of Bal with |0 > as F=2 state) into je > state (a hyperfine 
sublevel F=4 of the N=l rotational state) using a microwave 
<t + pulse. The second microwave pulse, coupling the J = 1 
and J — 2 states serves as a dressing field, preventing from 
non-radiative energy exchange between two interacting as well 
as other molecules, which are in the ground N = state. 



|00) - |00>, 
|01) - |01), 
110} - |io} , 

111) - *lll>. 

The resulting gate gives only " half of the true phase 
gate, so we need to perform the operation two times to 
obtain the phase gate, in which |11) changes sign. 

Electric fields required to mix the J = and J = 1 
rotational states of the ground state in order for the 
J = state to acquire a dipole moment will produce 
an interaction term /j,£/Ti ~ 7 • 10 10 s _1 , exceeding the 
dipole-dipole interaction strength u between two neigh- 
boring molecules, calculated above, by five orders of mag- 
nitude. It means that the electrostatic interaction contri- 
bution to the phase given in Eqs.©, (fTD]) would exceed 
the dipole-dipole interaction contribution from more than 
10 5 molecules in the lattice, and $ Cjt , <J> t « ii£T 27r /'h. In 
order to keep the phase error below %1 level, the elec- 
tric field magnitude will have to be controlled with better 
than 5 /xV/cm precision for = 100 /is. 



V. ROTATIONAL SCHEME BASED ON Bal 

To implement the rotational scheme we can use alkaline 
earth monohalides each having X 2 S + ground state with 
a large dipole moment and significant hyperfine split- 
tings. Specifically we consider 138 Ba 127 I, 138 Ba 79 Br and 
88 Sr 127 I, 88 Sr 79 Br with permanent ground state dipole 
moments [i rs 5.5 — 6 D [39(. 

A. Hyperfine structure of Bal 

Rotational, fine and hyperfine structure of the ground 
X 2 S+ state of 138 Ba 127 I can be calculated from the 
Hamiltonian 

H = H rot + H sr + Hhfs, (11) 

including the rotational H rot — B e N , spin-rotational 
H sr = J sr NS and hyperfine Hhf s — bpIS + d z >S z i — 
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eQq 3/ 2 - 1(1 + 1) /4J(7- 1) terms [37J. The coupling 

constants are j sr = 75.85 MHz, b F = 93.117 MHz, c = 
52.17 MHz and eQq = 33.62 MHz [13 ■ The spin-rotation 
and hyperfine terms are comparable, and the coupling of 
rotational, electron and nuclear spin angular momenta 
in the lowest rotational states is intermediate between 
the Hund's bp j and bps cases |38j . The energies of the 
hyp erfine levels were calculated using bpj coupling case 
[l7| . The hyperfine structure is due to the nuclear spin of 
127 I 1=5/2 (100% abundance), since the most abundant 
138 Ba isotope (71%) has zero nuclear spin. The scheme 
of the lowest rotational states of Bal including hyperfine 
structure is given in Fig(6] 



B. Phase gate operation time 

The qubit |0) and |1) states can be encoded into the 
F = 2 and F = 3 hyperfine sublevels of the ground rovi- 
brational state. The state |e) can be a hyperfine sub- 
level of either the N = 1 or N = 2 rotational state. 
On the one hand, it is preferable to choose the N — 2 
state since there is no electric-dipole allowed transition to 
the ground N — state, meaning that non-radiative en- 
ergy exchange processes between pairs of molecules due 
to electric dipole-dipole coupling, one in N = and the 
other in TV = 2 states are prohibited. On the other hand, 
the molecules can be excited to the N = 2 state only via 
two-photon processes, always resulting in a tt phase shift 
of the wavefunction (not the desired it/2 phase shift), 
which does not allow one to realize the phase gate. The 
solution can be to use a superposition of N = 1 and 
N = 2 levels, produced by a microwave field resonantly 
coupling these states, as shown in Fig[7^,. 

The resulting dressed states |+) and |— } are shifted 
from the unperturbed N — 1 level position by A Cit = 
±^ c ,t coupi , which is the effective Rabi frequency of cou- 
pling fields. The non-radiative energy exchange between 
two molecules has a strength Vind-dip ~ fi?/r 3 h, where 
/j, is the dipole moment of the corresponding transition, 
r is the distance between the molecules. For rotational 
transitions we consider the dipole moment is given by 
the permanent dipole moment of the ground electronic 
state. If |A Cjt | ^> Vi n d-dip the energy exchange be- 
tween the target and the control molecule as well as 
with other molecules in N — states is strongly sup- 
pressed, since the two molecules are no longer in reso- 
nance. Here we should note that for transitions between 
rotational levels of the same electronic and vibrational 
state the transition dipole moment is given by the per- 
manent dipole moment of the ground electronic state, 
so that the natural linewidth of a rotational level can 
be roughly calculated as A = 4^ 2 (27r) 3 /3A 3 rlu ,?i, where 
X mw is the wavelength of the N — > N — 1 transition. 
Since \ mw is larger or comparable to a typical distance 
r between molecules in an optical lattice or in electro- 
static traps, A < Vi n d~dip and the exchange processes 
are strong. For the case of CO and LiCs analysed in 



previous sections, r > A/2 meaning that A > Vi n d-dip 
and the exchange processes are slower than the spon- 
taneous decay and can be neglected. In the case of Bal 
the \N = 1,J= 3/2, F = 4} and \N = 2, J = 5/2, F = 5) 
states can be mixed to form |+) and |— ) states in the con- 
trol and target qubits. FigJTJa illustrates the way qubits 
are encoded in hyperfine sublevels of the N = state aa 
well as the coupling of the N = 1 and N = 2 states to 
form the \e) state. The phase gate in this case is direct 
since the ground N = state does not possess a dipole 
moment. The gate can be implemented in the same way 
as in the CO case. First, a tt pulse excites the control 
qubit into |e) = |+) state. Next a 2ir pulse, resonant with 
the unshifted |1) — > |+) = |e) transition is applied to the 
target molecule. If the control qubit is in the |e) state, 
the doubly excited |ee) state will be shifted in energy by 
the Vdip — /i 2 /r 3 h and the 2ir pulse will be out of reso- 
nance, producing a small phase shift of the wavefunction 
significantly less than tt. Finally, the control molecule is 
returned back to |1) by a second it pulse. 

As was already discussed above, Raman pulses used to 
excite a molecule from |1) to |e) always bring about a 7r 
shift of the wavefunction, which does not allow one to 
realize the phase gate. In turn, it means that the rota- 
tional scheme can be implemented only using one-photon 
microwave pulses, directly driving rotational transitions, 
making this scheme suitable for the architecture involv- 
ing electrostatic traps and coupled microwave resonators 
14j . Let us estimate the time required to implement 
the gate. As was shown the splitting of the dressed 
states |+) and |— ) has to exceed the exchange interac- 
tion strength Vind-dip — [i 2 /h 2 rfi (for coupled resonators 
architecture). Taking /j, « 6 D, h = 0.1 fim, r = 10 
/xm we have Vind-dip ~ 3.6 ■ 10 5 s . So we can take 
&>c,cou P i ~ 3 ■ 10 7 s -1 and tt t ,cou P i ~ 3 ■ 10 6 s~\ corre- 
sponding to intensities of the coupling MW fields I c ~ 6 
mW/cm 2 and It ~ 60 /iW/cm 2 , respectively. The Rabi 
frequency of the it pulses used to manipulate the control 
molecule has to be less than A c to be able to couple the 
ground state to only one of the dressed states. We can 
take Cl n = 3 • 10 5 s -1 , resulting in the duration of the 
pulse = ir/Qn = 10 fis. The target qubit 2tt pulse 

ind—dip-, SO we 



Rabi frequency has to be "C Vdip ~ Vj 



can take = 2 • 10 s . The resulting duration of the 
2tt pulse is T 2 tt = 2tt/Q 2 tt ~ 300 /its. So the total phase 
gate time is T gate = 320 [is. This time is significantly 
smaller than the natural decay time of rotational states 
of the ground electronic and vibrational state of a polar 
molecule (typical numbers are 10 5 s). 

The state of the qubit can be read out using a mi- 
crowave field tuned far from resonance with the qubit 
transition [l4j. After the field interacts with the qubit 
transition, it accumulates a phase, which depends on the 
state of the qubit, allowing one to read out the qubit 
state from the transmission or reflection spectra of the 
field. 

The ground rovibrational state N = of Bal is strong- 
field seeking, and there are proposals as to how to de- 
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FIG. 7: a) Proposed scheme of suppression of non-radiative 
energy exchange between a pair of molecules; b) Hyperfine 
sublevels of the ground rotational state N = used to encode 
qubits in the control and target molecules. 



time of a molecule in an optical lattice and an electro- 
static trap. Lifetime of a nuclear spin state of a single 
molecule, isolated from the environment, can be as long 
as years. Lifetimes of ultracold molecules in a far-detuned 
optical lattice of « 1 s have been obtained by minimizing 
scattering of lattice photons [Z(| • 

Since all three schemes rely on the dipole-blockade 
mechanism, it is not critical (although it is desirable) to 
cool molecules to the motional ground state of a lattice or 
an electrostatic trap. It is experimentally possible to load 
BEC atoms, from which molecules can be formed using 
photoassociation or Feshbach resonance techniques, into 
the ground state of a lattice by switching it on adiabati- 
cally [il[ . In an electrostatic trap sideband cooling using 
microwave resonator enhanced spontaneous emission was 
proposed for efficient cooling of molecules to the ground 
state of motion [14j . With overall coherence time of the 
order of 1 s, the maximal number of operations, which 
can be performed with the direct, inverted and rotational 
gates, is ~ 8 ■ 10 3 , 6 • 10 3 and 3 • 10 3 , respectively. 



sign an electrostatic trap of this type. For a pair of 
molecules, each trapped above a small conductor, with 
the conductors connected by a superconducting wire, 
the dipole-dipole interaction is modified and has a form 
Vdip ~ [J, 2 /h 2 r. Actually, this is strictly true in the 
case when the temporal dynamics of interacting dipoles 
is slower than the characteristic frequency of the wire, 
scaling as lo = nirv/r, where n is an integer, and v is the 
transmission velocity of the wire. Taking r ~ 10 /xm, and 
v ~ c, where c is the speed of light, we have the lowest 
frequency of the wire tu m i n ~ 10 13 s -1 . It corresponds to 
characteristic times ~ 100 fs, significantly smaller than 
any processes in our system. It means that the charge dis- 
tribution in the wire will adiabatically follow the dynam- 
ics of the dipoles, and the coupling between the dipoles 
can be considered as static. The mixing between N = 1 
and N = 2 states can be realized in a microwave cavity 
as was proposed in Ref . fl4j , the coupling strength can be 
varied by tuning the frequency of the N = 1 «-> N = 2 
transition in and out of resonance with the cavity mode 
by the trap DC electric field. Microwave ir and 2-7T pulses 
can be realized using classical microwave fields, resonant 
to the |1) — ► \e) transition of a selected molecule. 

VI. DECOHERENCE MECHANISMS 

In the proposed schemes, qubits are stored in hyperfine 
sublevels of a ground rovibrational electronic state, which 
makes them insensitive (at least in first order) to local 
fluctuations of DC and AC electric fields. The sublevels 
do feel fluctuations in a magnetic field, which should be 
minimized on a time scale of a second, relevant to a life- 



VII. CONCLUSION 

Analysis of quantum computation proposals with polar 
molecules allows one to define a set of requirements for 
a molecular system optimizing its performance in terms 
of quantum gate time, number of operations, coherence 
time, gate error and experimental simplicity. A scheme 
studied in detail in the present work utilizes switchable 
dipole-dipole interaction between molecules, which can 
be used to implement a universal two-qubit gate such as 
a phase gate. It was demonstrated that proposed schemes 
of the phase gate realization using polar molecules sat- 
isfy the set of requirements presented in the introduction, 
and are experimentally feasible with currently existing 
technologies. For the direct phase gate the CO and NF 
molecules are suitable candidates, and we are looking for 
other molecules with the same parameters, while for the 
inverted and rotational schemes two classes of molecules 
are suggested, namely, alkali dimers and alkaline-earth 
monohalides, respectively. Optical lattice and architec- 
tures using electrostatic traps and coupled microwave res- 
onators combined with the dipole-blockade effect can be 
used to implement the phase gate with high fidelity. De- 
coherence present in both architectures limits the qubit 
lifetime to ~ 1 s, allowing several thousand operations to 
be performed. 
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